Journal  of  Coastal  Research 


SI 

59 

1-6 

West  Palm  Beach,  Florida 

2011 

The  Coastal  Modeling  System  Flow  Model  (CMS-Flow):  Past  and 
Present 


Christopher  W.  Reed  ,  Mitchell  E.  Brown*,  Alejandro  Sanchez*,  Weiming  \Vus,  and  Adele  M.  Buttolph 

It  id  c  CnmArotmn  §National  Center  for  Computational  Hydroscience  and 

Engineering, 

University  of  Mississippi, 


TURS  Corporation 
1625  Summit  Lake  Drive 
Tallahassee,  FL  32317,  USA 


Chris  _reed@urscorp  .com 


University,  MS,  28677,  USA 


JU.S.  Army  Engineer  Research  &  Development  Center,  ^Coastal  Analysis  LLC, 

3909Halls  Ferry  Road  29  Valley  Heights  Drive, 

Vicksburg,  MS  39180-6199,  USA  Enterprise,  OR  97828,  USA 


www.  cerf-jcr.  org 


ABSTRACT 

REED,  C.W.;  BROWN,  M.E.;  SANCHEZ,  A.;  WU,  W.,  and  BUTTOLPH,  A.M.,  2011.  The  Coastal  Modeling 
System  Flow  Model  (CMS-Flow):  Past  and  Present.  In:  Roberts,  T.M.,  Rosati,  J.D.,  and  Wang,  P.  (eds.), 
Proceedings,  Symposium  to  Honor  Dr.  Nicholas  C.  Kraus,  Journal  of  Coastal  Research,  Special  Issue,  No.  59,  pp.  1- 
6.  West  Palm  Beach  (Florida),  ISSN  0749-0208. 

CMS-Flow  is  a  coupled  time-dependent  circulation,  sediment  transport  and  morphodynamic  model  based  on  the 
numerical  solution  of  the  mass,  momentum  and  transport  equations  on  a  Cartesian  (quad-tree)  grid  network  with  both 
explicit  and  implicit  solvers.  It  has  been  developed  and  is  currently  supported  under  the  Coastal  Inlets  Research 
Program  (CIRP)  conducted  at  the  U.S.  Army  Engineer  Research  and  Development  Center  (ERDC),  Coastal  and 
Hydraulics  Laboratory  (CHL).  The  model’s  primary  function  is  to  support  multi-disciplinary  research  teams  and 
conduct  practical  projects  at  coastal  inlets.  CMS-Flow  has  been  designed  with  a  relatively  simple  code  structure 
which  allows  for  rapid  development  and  inclusion  of  new  sediment  transport  algorithms,  while  always  being 
accessible  to  the  general  modeling  community,  including  both  USACE  and  commercial  users.  Today,  CMS-Flow  is 
an  integral  component  of  the  CIRP,  providing  technology  for  simulating  hydrodynamics,  waves,  sediment  transport 
and  morphology  for  short  and  long  timeframes  in  coastal  inlets,  adjacent  beaches,  navigation  channels  and  bays. 

ADDITIONAL  INDEX  WORDS:  Coastal  hydrodynamics,  sediment  transport,  morphodynamic  modeling. 


INTRODUCTION 

CMS-Flow  is  a  coupled  time-dependent  circulation,  sediment 
transport  and  morphodynamic  model  that  can  be  configured  for 
both  2D  and  3D  simulations.  The  model  is  based  on  the 
numerical  solution  of  the  mass,  momentum  and  transport 
equations  on  a  Cartesian  grid  network  with  both  explicit  and 
implicit  solution  schemes  available.  The  circulation  model 
component  supports  spatially  variable  and  time  dependent  wind, 
tide,  atmospheric  pressure,  flow  and  wave  radiation  stress 
forcing,  salinity  transport  and  includes  special  features  such  as 
flooding  and  drying,  spatially  variable  bottom  friction, 
representation  of  flow  control  structures,  and  advanced 
boundary  condition  treatments.  Sediment  transport  is  forced  by 
the  circulation  model  derived  currents  and  waves  derived  from 
the  coupled  CMS-Wave  model  (Lin  et  al.,  this  issue).  Various 
bedload  and  suspended  load  models  such  as  Lund-CIRP  and 
Watanabe  formulations  are  available  via  user  options.  Recently, 
a  non-equilibrium  transport  formulation  was  included  in  the 
model  (Sanchez  et  al. ,  this  issue)  and  sediment  avalanching, 
hard-bottom  and  breaching  algorithms  were  added.  The 
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hydrodynamic  and  sediment  transport  models  are  coupled  via  a 
morphodynamic  module  that  calculates  time-dependent  changes 
in  the  bed  elevation  due  to  the  simulated  sediment  transport.  The 
hydrodynamic  grid  is  then  updated  to  reflect  the  depth  changes, 
which  then  completes  the  coupling. 

CMS-Flow  has  been  developed  and  is  currently  supported 
under  the  Coastal  Inlets  Research  Program  (CIRP)  conducted  at 
the  U.S.  Army  Engineer  Research  and  Development  Center 
(ERDC),  Coastal  and  Hydraulics  Laboratory  (CHL).  The 
model’s  primary  function  is  to  support  multi-disciplinary 
research  teams  and  conduct  practical  projects  at  coastal  inlets. 
CMS-Flow  has  been  applied  at  numerous  coastal  inlets  and 
examined  extensively  through  sensitivity  testing.  It  has  been 
designed  with  a  relatively  simple  code  structure  which  allows 
for  rapid  development  and  inclusion  of  new  sediment  transport 
technologies,  while  always  being  accessible  to  the  general 
modeling  community,  including  both  USACE  and  commercial 
users. 

CMS-Flow  is  supported  with  the  Surfacewater  Modeling 
System  (SMS)  software  package,  which  provides  a  graphical 
user  interface  for  model  setup,  running  and  post-processing. 
Details  of  the  theory  and  numerical  implementation  can  be 
found  in  the  CMS-Flow  technical  reports  (Milletello  et  al., 
2004;  Buttolph  et  al.,  2006)  and  a  description  of  the  graphical 
interface  can  be  found  in  the  SMS  literature  (Zundel,  2000). 
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BACKGROUND 

CMS-Flow  began  in  1994  as  a  ID  hydrodynamic  model, 
developed  by  Adele  Militello  as  a  class  project  taught  by  Dr. 
Don  Resio  at  Florida  Institute  of  Technology  and  was  called 
M2D.  With  guidance  by  her  advisor,  Dr.  Gary  Zarillo,  and 
continued  support  from  Dr.  Resio,  the  model  was  extended  to 
two  dimensions.  The  numerical  solution  was  explicit  and 
employed  a  Cartesian  grid,  both  relatively  simple  numerical 
approaches  which  today  still  represent  the  underlying 
philosophy  behind  CMS-Flow:  to  provide  a  modeling  system 
that  would  be  easy  for  users  to  set  up  and  apply,  code  that  was 
readily  understandable  for  developers,  code  that  was  easy  to 
work  within  and  modify,  and  excellent  performance  on 
numerical  integrity  tests.  CMS-Flow  is  designed  to  provide  an 
efficient  means  to  develop  and  test  new  transport  algorithms  and 
allow  for  validated  technologies  to  be  immediately  available  to 
the  user  community. 

While  studying  at  Florida  Institute  of  Technology,  Militello 
took  a  position  at  Texas  A&M  University  at  Corpus  Christi  to 
work  with  Dr.  Nick  Kraus.  The  first  application  of  the  M2D 
model  was  to  the  Lower  Laguna  Madre,  Texas.  However,  no 
interface  for  grid  development  was  available,  so  the  numerical 
grid  was  completely  hand-edited  after  being  sketched  out  on 
about  8  feet  of  paper,  and  took  about  three  weeks  to  complete. 

In  the  latter  part  of  the  1990's,  the  CIRP  program  was 
evaluating  approaches  to  coastal  hydrodynamic  modeling.  When 
Dr.  Kraus  became  the  director  of  CIRP,  he  immediately 
recognized  the  benefits  of  the  M2D  model  attributes  and 
adopted  it  into  the  CIRP  program  as  the  primary  hub  for 
advancing  sediment  transport  prediction  technology.  The  M2D 
model  was  eventually  renamed  CMS-Flow  for  Coastal  Modeling 
System  to  better  reflect  its  emerging  role  in  coastal  modeling, 
and  an  interface  within  the  SMS  package  was  created.  Model 
enhancements  have  steadily  been  added  to  CMS-Flow  as  well  as 
to  its  companion  wave  model  CMS- Wave,  providing  the 
modeling  community  with  a  robust,  stable  and  technologically 
advanced  modeling  system  that  can  be  used  efficiently  by  both 
expert  modelers  as  well  as  those  with  expertise  in  sediments  and 
related  fields,  but  not  necessarily  computational  fluid 
mechanics.  Today,  the  CMS-Flow  software  is  an  integral 
component  of  the  CIRP,  providing  technology  for  simulating 
hydrodynamics,  waves,  sediment  transport  and  morphology  for 
short  and  long  timeframes  in  coastal  inlets,  adjacent  beaches, 
navigation  channels  and  bays.  The  CIRP  team  comprises 
numerous  contributors  including  scientists  and  engineers  with 
backgrounds  in  coastal  sciences,  geology,  sediment  transport, 
numerical  methods  and  software  development.  CMS-Flow  has 
been  used  in  over  80  projects  within  the  U.S.  and  more 
internationally.  The  project  locations  of  both  US  ACE  and 
commercial  applications  are  shown  in  Figure  1. 

MODEL  DESCRIPTION 

The  CMS-Flow  2D  model  equations  are  based  on  the  shallow 
water  equations.  The  2D,  depth-integrated  continuity  and 
momentum  equations  solved  by  CMS-flow  are 
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dt  dx  dy 
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where  the  subscripts  x  and  y  indicate  the  x  and  y  components 
and: 


h  =  still-water  depth 

r|  =  water- surface  elevation 
t  =  time 

u  =  depth-averaged  velocity,  x- 
direction 

v  =  depth-averaged  velocity,  y- 
direction 

xb=  bottom  stress 


g  =  acceleration  due  to 
gravity 

f  =  Coriolis  parameter 
q  =  flow  per  unit  width 
D  =  diffusion  coefficient 

xw  =  surface  stress 

xs  =  wave  stress 


A  variety  of  sediment  transport  formulations  are  available. 
For  instance,  the  Advection-Diffusion  (AD)  scheme  uses  the 
Lund-CIRP  bed  load  and  a  suspended  load  transport  equation. 
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The  bed  load  formulation  was  developed  by  Camenen  and 
Larson  (2005): 


Q be 


J(s-l)gdl 


Qc 


yCW' 


(4) 


where 

qbc  =  bed  load 


d50=  median  grain  size 


0. 


=  mean  wave-current  an,b  =  model  coefficients 


Shields  parameter  stress 

0C  =  current  Shields  parameter  s  =  specific  gravity  of 

sediment 

9CW  =  maximum  wave-current  0cr  =  critical  Shields  parameter 
Shields  parameter 


and  the  AD  equation  is: 


d(Cd)  d(Cqx )  d(Cqy) 
dt  +  dx  +  dy 


where 

C  =  depth-averaged  sediment  K  =  sediment  diffusion 
concentration  coefficient 

d  =  total  water  depth  D  =  sediment  deposition 

rate 


P  =  sediment  pick-up  rate 

The  sediment  continuity  equation  is  solved  to  compute  depth 
change 


dh 

dt 


i  r  9Qx 

pwg  \l-pj  V  dx 


(6) 


where  Qx  and  Qy  are  bed  load  in  the  x  and  y  directions,  p  is  the 
sediment  porosity  and  h  is  the  bottom  elevation. 

This  brief  description  of  the  governing  equations  provides  the 
general  mathematical  basis  for  the  CMS  model.  The  detailed 
description  of  the  governing  equations,  parameter  definitions, 
numerical  solution  schemes  and  additional  algorithms  are 
available  in  the  referenced  report.  The  formulations  for  the  3D 
module  are  available  in  the  draft  model  report  (Reed  et  al ., 
2007). 


RECENT  MODEL  ENHANCEMENTS 

Three  key  components  inherent  in  choosing  and  developing  a 
model  in  practice  are: 

•  Capabilities 

•  Reliability 

•  Computational  Speed 


In  any  coastal  modeling  application  it  is  critical  to  use  a 
model  that  can  represent  the  important  and  controlling  physical 
processes  and  controls  on  the  system.  These  can  range  from 
fundamental  processes  such  as  cohesion  to  the  ability  to 
represent  spatially  variable  wind  effects.  It  is  also  important  that 
the  model  is  reliable,  and  can  be  setup  and  applied  with  the 
confidence  that  a  stable  simulation  can  be  obtained  without 
spurious  results  or  the  need  for  tedious  and  continual  grid 
modifications.  Finally,  despite  the  processing  speeds  of  modem 
computers,  the  desire  to  include  more  refined  models  with  more 
detailed  physics  has,  and  always  will,  push  the  limits  of 
available  computational  resources. 

Enhancements  to  the  CMS-Flow  model  are  made  regularly  to 
improve  the  model  capabilities,  reliability  and  computational 
efficiency.  Typically  in  model  development,  the  time  required 
for  a  new  model  feature  to  be  supported  by  the  GUI  (i.e.  SMS) 
so  that  the  general  users  can  access  the  features  may  take  six 
months  to  a  year.  In  order  to  reduce  the  time  for  new  features  to 
become  available,  the  advanced  card  option  has  been 
implemented  in  the  SMS  GUI.  The  advanced  card  option  allows 
the  user  to  enter  special  commands  to  the  CMS-Flow  model  in  a 
special  SMS  menu  and  these  commands  will  be  read  by  the 
CMS-Flow  code  and  control  the  model  response.  The  menu 
structure  is  static,  but  will  accommodate  new  commands  as  they 
are  added  to  the  model  code.  These  features  allow  the  code 
developers  to  add  new  or  experimental  features  without 
requiring  a  corresponding  update  to  the  SMS.  The  advanced 
cards  have  been  used  to  add  additional  friction  algorithms, 
diffusion  algorithms,  avalanching  schemes,  flow  structures, 
multiple  grain  size  simulations,  and  multiple  processor 
capabilities.  Eventually,  the  capabilities  added  with  the 
advanced  cards  are  supported  by  SMS.  An  example  of  recent 
capabilities  added  to  the  CMS-Flow  using  the  advanced  cards 
are  a  variety  of  flow  structures,  including  weirs,  culverts, 
pilings,  tidal  gates  and  porous  rubble  mounds. 

An  example  of  a  model  enhancement  developed  to  improve 
the  model  reliability  is  the  wave  adjusted  boundary  condition 
(Reed,  2005).  In  the  nearshore,  the  wave-induced  longshore 
current  and  setup  are  dominant  hydrodynamic  processes,  and 
lateral  boundary  conditions  representing  tide  and  oceanic 
forcing  will  not  include  corresponding  surface-wave 
contributions.  Without  proper  boundary  conditions,  significant 
gradients  in  current  and  water  level  can  occur  that  contaminate 
the  solution  in  the  internal  domain.  This  situation  is  depicted  in 
Figure  2.  A  standard  strategy  in  such  situations  is  to  place  the 
boundaries  far  from  the  site  of  interest,  but  this  strategy  greatly 
increases  computer  resource  requirements  and  run  time,  and  it 
may  not  be  appropriate  for  long-term  simulations.  A  wave- 
adjusted  boundary  condition  that  accounts  for  wave-induced 
water  level  and  current  acting  in  combination  with  tidal  forcing 
has  been  developed  and  implemented  in  CMS-Flow.  The 
approach  consists  of  solving  the  one-dimensional  (ID) 
continuity  and  momentum  equations  containing  the  radiation- 
stress  gradients  along  the  axis  normal  to  the  boundary.  Wave- 
induced  setup  along  this  axis  is  added  to  the  prescribed  water- 
elevation  boundary  forcing  values  to  provide  consistent  water 
elevation  values  at  the  boundary  and  in  the  grid  interior,  which 
subsequently  eliminates  anomalous  velocities.  The  effect  of  the 
wave  adjusted  boundary  conditions  have  been  applied  to  the 
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Figure  2.  Velocity  calculated  with  standard  water  surface  elevation 
boundary  condition. 


Figure  3.  Velocity  calculated  with  wave-adjusted  boundary  condition. 


simulation  depicted  in  Figure  2  and  the  results  are  shown  in 
Figure  3.  The  anomalous  velocities  at  the  boundary  have  been 
eliminated  by  the  wave  adjusted  boundary  conditions. 

A  number  of  enhancements  to  CMS-Flow  have  been  made  to 
increase  the  computational  efficiency.  The  most  direct  approach 
is  the  application  of  OpenMP  for  use  on  multiple  processors 
with  shared  memory.  This  enhancement  consisted  of 
restructuring  portions  of  the  code  to  allow  for  OpenMP 
processing  and  adding  the  OpenMP  directives. 

A  second  enhancement  is  the  additional  of  an  implicit  solver. 
Often  in  applications  of  the  CMS-Flow  model  with  the  explicit 
solver,  the  time  step  is  limited  by  numerical  stability  criteria 
rather  by  the  temporal  scale  of  the  physical  processes.  The 
implicit  solver  eliminates  the  numerical  restriction,  and  although 
it  requires  more  computational  time  per  time  step,  it  can  use  a 
significantly  larger  time  step  and  decrease  simulation  CPU  time 


by  an  order  of  magnitude  for  many  applications.  However, 
adding  new  capabilities  within  the  implicit  solver  requires  more 
development  time  and  some  processes,  such  as  wetting  and 
drying,  require  smaller  time  steps  making  the  explicit  solver  the 
more  efficient  option.  Therefore,  the  CIRP  program  will 
maintain  both  an  explicit  and  implicit  solver,  and  the  user  can 
select  the  solver  most  suitable  for  the  specific  application. 

The  third  enhancement  to  improve  computational  efficiency  is 
the  incorporation  of  the  quad-tree  grid  structure  (Wu  et  al. ,  this 
issue).  The  original  Cartesian  grid  structure  has  provided  for  a 
simple  solution  algorithm  and  a  solid  basis  for  rapid  code 
development.  However,  the  grid  structure  can  be  cumbersome 
when  trying  to  resolve  complex  geometries  over  large  areas.  The 
structure  requires  a  large  number  of  grid  cells  and  for  the 
explicit  solver  will  also  yield  additional  time  step  restrictions 
due  to  the  inherent  extension  of  fine  grid  cells  to  deep  water 
locations.  The  quad-tree  grid  structure  is  a  non-conforming  grid 
with  a  maximum  of  two-to-one  cell  splitting  along  any  cell  face. 
The  quad-tree  structure  will  reduce  the  total  number  of  grid  cells 
needed  to  resolve  a  model  domain,  thus  benefiting  both  the 
explicit  and  implicit  CMS  solvers.  An  additional  benefit  for  the 
explicit  solver  is  the  increase  in  time  steps  allowed  since  small 
grid  cells  will  not  be  extended  into  deep  water.  An  example  of 
the  benefits  of  using  the  quad-tree  are  shown  in  Figures  4,  5  and 
6  show  an  example  of  the  benefits  of  using  the  quad-tree  grid 
structure.  An  example  model  domain  in  Cape  Sable,  Florida  is 
shown  in  Figure  4,  and  contains  a  series  of  small  canals 
embedded  in  an  intertidal  zone  adjacent  to  a  lake.  In  order  to 
fully  resolve  the  canals,  four  cells  within  the  canal  cross-section 
of  the  grid  shown  in  Figure  5  were  required.  The  grid  orientation 
was  adjusted  to  minimize  the  number  of  grid  cells.  For  this  grid 
a  total  of  57,000  cells  were  required.  The  time  step,  which  is 
controlled  by  the  shallow  water  wave  speed  in  the  deeper  water 
areas,  was  limited  by  the  small  cells  extending  to  the  west,  and 
was  limited  to  0.15  seconds.  A  quad-tree  grid  for  the  same 
domain  is  shown  in  Figure  6  and  contains  the  same  small  grid 
cell  size  within  the  canals.  This  grid  contains  only  2000  cells, 
and  the  allowable  time  step  is  0.4  seconds.  This  combined 
computational  speed  increase  for  an  explicit  solver  due  to  both 
the  reduced  cell  count  and  larger  time  step  is  a  factor  of  over  25. 


Figure  4.  Cape  Sable  Florida  Canal  System. 


Journal  of  Coastal  Research,  Special  Issue  No.  59,  2011 


Coastal  Modeling  System  Flow  Model  (CMS-Flow) 


5 


MODEL  VALIDATIONS 

Despite  the  push  for  rapid  development  and  application  of 
new  technology,  each  component  of  the  CMS-Flow  model  is 
rigorously  tested  prior  to  being  released  for  general  use.  The 
testing  includes  both  model  calibration  and  validation  using 
analytical,  laboratory  and  field  data  and  testing  on  a  set  of 
standard  bench  test  cases.  This  second  set  of  tests  is  required  to 
assure  that  the  new  modules  function  properly  within  all  existing 
forcing  and  boundary  conditions  and  that  their  implementation 
has  not  compromised  existing  modules  and  functionality.  The 
classic  calibrations  and  validations  are  used  to  set  model 
parameter  values,  or  ranges  for  those  parameter  values,  and 
assesses  the  model  accuracy  and  reliability. 


As  an  example,  Figure  7  shows  the  validation  of  the  CMS- 
Flow  model  response  to  wave  forcing,  calculations  were 
compared  to  laboratory  data  collected  by  Visser  (1982).  The 
CMS-Flow  model  was  coupled  with  the  CMS- Wave  model 
using  the  SMS  interface  to  develop  these  simulations.  This 
calibration  was  used  to  validate  the  surf-zone  diffusion  module. 

As  a  second  example,  Figure  8  shows  the  comparison  of  the 
model  simulation  of  cross-channel  transport  and  morphology 
changes  compared  to  the  data  of  Galappatti  and  Vreugdenhil 
(1985).  The  test  case  represents  a  51  cm/s  cross-channel  flow 
over  a  39  cm  deep  channel  with  an  initial  side  slope  of  1:10  and 
0.16  mm  sand  particles.  These  comparisons  were  used  to 
determine  the  adaptation  length  scale  in  the  non-equilibrium 
sediment  transport  formulation. 


Figure  5.  Cartesian  Grid  covering  canals  and  opening  to  lake. 


Figure  6.  Quad-tree  grid  covering  canals  and  opening  to  lake. 


Figure  7.  Comparison  of  CMS-Flow  to  data  of  Visser  (1985). 


Figure  8  Calibration  of  the  Non-Equilibrium  Sediment  Transport 
Algorithm. 


SUMMARY 

The  CMS-Flow  model  has  become  an  integral  component  of 
the  CIRP.  When  coupled  with  the  companion  wave  model  CMS- 
Wave  (Lin  et  al.,  this  issue)  the  CMS  model  provides  a  stable 
and  robust  computational-based  analysis  tool  for  simulating 
complex  sediment  transport  and  morphologic  response  of  inlets, 
beaches  and  navigation  channels  to  forcing  by  tides,  winds,  river 
flows  and  waves. 
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